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We introduce a new sum-rule for large- N c QCD which relates the den- 
| 7 | \ sity of heavy quarkonium states, the state-averaged square of the wave- 

£f~) \ function at the origin, and the heavy quark current- current correlator. 

. Focusing on the region of energy just above perturbative threshold, we 

calculate the correlator by incorporating arbitrarily high orders in the 
CN ■ QCD coupling a s . We use the sum-rule to determine the bottomo- 

O ' nium potential using experimentally measured s-wave leptonic widths 

^ | and compare the result with the potential obtained by direct calcula- 

tion from the measured s-wave spectrum. We discuss the utility of the 



sum-rule method for accurate determination of the confining potential. 



^ 1 Introduction 

> : 

^ ■ In the heavy quarkonium system, we can study quark confinement without the added com- 
plexity of relativistic processes. As the mass of the quark becomes large, the dynamics of 
the quarkonium system becomes non-relativistic and can be described by the Schrodinger 
equation with a static potential. The short distance behavior of this potential can be ap- 
proximated by single gluon exchange and therefore has the familiar Coulombic form. In this 
paper, we consider the remaining long distance part of the potential, which we refer to as 
the confining potential. In the course of our analysis we answer the following questions: (i) 
Can the confining potential in some region of r be determined from the spectrum by direct 
calculation, i.e., without parameter fits? (ii) Can the confining potential also be determined 
from other data, such as the s-wave leptonic widths T e + e -7 (iii) How do the results of (i) 
and (ii) compare, and can we accommodate both data sets with one potential? 

1 Research supported in part by the National Science Foundation under Grant #PHY-9218167. 
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To answer these questions we make use of the large- N c expansion, where N c is the number 
colors. We derive a sum-rule that relates the imaginary part of the heavy quark current- 
current correlator with the product of p, the density of s-wave states, and |^ e (0)| , the 
state-averaged square of the wavefunction at the origin. The derivation involves calculating 
the imaginary part of an integral of the heavy quark current-current correlator in momentum 
space. We do this calculation in two different ways, once using perturbative methods and 
then again using spectral sums. The sum-rule follows from equating the two calculations. 

Sum-rules for the current-current correlator for light quark systems have a long his- 
tory [|J. Our approach is new in that we restrict our attention to heavy quark currents 
and focus on the energy region just above the perturbative qq threshold. As it turns out 
this energy region has some interesting characteristics. One feature is the appearance of 
threshold singularities in the perturbative calculation. To account for these singularities, 
our calculation will include contributions from arbitrarily large orders in a s , unlike more 
conventional treatments. Another feature is that the effects of confinement can be described 
using non-relativistic bound state methods. As we attempt to illustrate, application of the 
sum-rule provides a useful complement to non-relativistic methods. In particular, we discuss 
its use in the determination of the confining potential for the bottomonium system. 

The organization of the paper is as follows. In the sections 2 and 3, we derive the sum- 
rule. In section 2, we discuss the perturbative calculation of the current- current correlator, 
and we note that it is straightforward to include information from arbitrarily high orders in 
a s by using the energy spectrum of the Coulomb bound states. In section 3, we calculate the 
current-current correlator as a spectral sum and relate it to the properties of the physical qq 
states in the large- iV c limit. In section 4, we test the sum-rule for the bottomonium system. In 
the remaining sections we discuss various ways to determine the confining potential, discuss 
our various approximations and the physics we have left out, and answer the questions (i), 
(ii), and (iii). 



2 Perturbative calculation 

We consider quantum chromo dynamics with N c colors and one flavor of heavy quark, which 
we call q. Let m be the mass of the heavy quark and J M (x) be the corresponding electric 
current. We can write the renormalized current-current correlator in momentum-space as 

/ d 4 xe ik * (0|r[J"(a:)J*(0)]|0)= hm ^ - ^-^U(Vk^,e), (1) 

where e is the infinitesimal positive parameter used to regulate the singularities of the bare 
propagators. Because we are in the time-like momentum region k 2 > 0, k° > 0, we can move 
to the rest frame where k = and rewrite U(\/k 2 ~ } e) as U(k°, e). 
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Let z and Az be fixed positive numbers such that Az << z << 1. Let us define the 
integral 

P(m,z,Az)= \im J n((2 + C)m, e) d£, (2) 

(eC(z,Az) 

where C is the contour in the complex ^-plane shown in Figure 1|] The purpose of this 
contour integral will be explained in a moment when we discuss the N c — > oo limit. Our 
choice of energy values places us within the non-relativistic qq resonance spectrum, just 
above the perturbative threshold, k° = 2m. 

If the dimensionless quantities z and Az are held fixed, we expect P to have a well- 
defined asymptotic expansion in a s (m 2 ) asm-> oo, where a s (m 2 ) is the QCD coupling^ at 
scale m 2 . In the following we abbreviate a s (m 2 ) as a s . We are interested in the imaginary 
(absorptive) part of P, which we write as 

o 

771 

Im P(m, z, Az) = N c — Az-G(z,a s ). (3) 

The perturbative expansion of G will be discussed below. Let us now divide out the explicit 
N c in (|J) and consider 

N' 1 Im P(m, z, Az) (4) 

in the simultaneous limit iV c — > oo, m — > oo. We rescale the coupling constant in the 
standard way so that a s Cp has a well-defined limit as N c — > oo. We then have 



lim 

N c — >oo 



iVT 1 hn.P(m,z, Az) 



2 

m 



Az-G(z,a s ). (5) 

Z7T 



As N c —> oo, the resonances become stable and the corresponding pole singularities 
approach the real £-axis from below. As m — > oo, we expect the density of such pole 
singularities approaching the real £-axis to increaseP] and anticipate that the integrand in (0) 
will exist in the double limit N c — > oo, m —> oo only in the sense of a distribution. The 
contour integral defining P, however, does not have such problems. The pole singularities 
approach the real axis from one side and therefore do not pinch the contour. This is the 
motivation for our definition of P. 

For a s z 1//2 , we can compute G in (^) by expanding in powers of a s , i.e., the usual 
perturbative loop expansion. This expansion, however, is inadequate for z 1 / 2 < a s . As the 



1 To avoid the possibility of endpoint singularities we have pulled the endpoints of the contour slightly off 
the real axis and into the upper half plane, i.e, we have added a small positive imaginary part to z. This 
technical point is of little consequence, however, since the limit Imz — > + presents no problems. 

2 As we will discuss section 6, the precise scale at which a s is evaluated does not matter to the order we 
are working. The difference between a s {m 2 ) and a s (zm 2 ) is not important. 

3 This is discussed in the next section. 
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order of a s increases, the order of the z 1 I 2 singularity also increases |3[]. It is therefore more 
useful to expand G as 

oo 

G{z,a s ) = X> S CV)" • G n {a s C F z~ 1 / 2 ), (6) 

n=l 

where Cp = J ^ c jV " 1 . The leading term, G\, is given by the sum of all simple ladder diagrams, 
which in turn can be obtained from the bound-state spectrum for the Coulomb potential. 
G\ has been calculated explicitly in || (see also ||), with the result 

We will use this form in section 4. 



3 Spectral sum calculation 

Let us now repeat the calculation in (j5|), this time using the spectral representation of the 
current-current correlator. We consider the contributions of all single and multiparticle states 
to the correlator: 

(0\T[J tM (x)J u (0)]\0) =6(x )J2(0\J^(x) \n) (n\r (0)\0) (8) 

n 

+ 6{-x°) ]T (0| J v (0) \n) (n\ J"{x) |0) . 

n 

In the N c — > oo limit, the coupling of the current to multiparticle meson states is suppressed 
by explicit factors of N c 2 . The contribution of these states to the spectral sum in (H) is 
therefore suppressed by iV" 1 , and we ignore such states in the iV c — >• oo limit. 

Having reduced the correlator to a spectral sum over qq resonance states, we use this 
spectral sum to compute the imaginary part of N~ X P. In the N c — > oo limit it is not 
necessary to include all resonance states. Resonances which do not approach the contour C 
give an imaginary contribution that is suppressed by AT" 1 ^ Let us use the index r to label 
resonances which approach the contour as N c — > oo. For each resonance r, let M r be the 
mass, r r be the width, and f^g^u be the coupling with the current. We then have 



N' 1 Im P(m, z, Az) = ^ Im 



-jvt 1 l/rl a dc 



+ O(N^) ■ (9) 



^From here on, we suppress writing the 0(N C x ) error term. We can compute the integral in 
terms of the residues of the poles which approach C, 



The real contribution, however, is not suppressed. 
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iV-WK;,A 2 H£|gg. (10) 

Since these resonances are non-relativistic bound states, we can use the van-Royen-Weisskopf 
result, 

Ar- 1 |/ r | 2 = 4|^ r (0)| 2 M r , (11) 

where Vv(0) is the Schrodinger wavefunction of the resonance at the origin. Only s-wave 
states couple to the photon in this limit. Substituting ([11]) into QlOD , we get 

AT 1 Im P(m, z, Az) = — ElVv(O) | 2 - (12) 

Comparing fl5|) and (|l|), we find 

-£lVv(0)| 2 = ^Az.G(*,a s )- (13) 

The spectral sum on the left side of ( O ) is a discrete sum. It therefore has intrinsic fluctu- 
ations of size [An(m)]" 1 , where An(zm) is the number of s-wave states between (2 + z)m 
and (2 + 2 + Az)m. These fluctuations are not reproduced in the perturbative calculation at 
any finite order in a s . This is what we expect from an operator product analysis in which the 
fluctuations are associated with matrix elements of higher dimension operators, proportional 
to powers of Aqcd H- We conclude that [An(zm)]" 1 is smaller than any power of a s as 
m — > 00. 

Let p(zm) be the density of s-wave states at energy zm. We can rewrite 

IVV(O) I 2 - P(zm) \^(0)\ 2 Az (14) 

III r 

where (^^(O)! 2 is the average of |?/v(0)| 2 for resonances with energies near zm. We then 
have 

2 

27rp(zm)|e(0)| 2 = ^G(^,a s ), (15) 

or 

P (zm) = m2G ^ a ^\ . (16) 
4^|«(0)| 2 

At this point, the careful reader may wonder whether we have done something incon- 
sistent, because we have argued in this section that the physical qq bound states affect the 
imaginary part of the current- current correlator only if they approach the contour C, while 
in section 2, we saw that the perturbative calculation of the same imaginary part receives 
contributions which appear to be associated with the Coulomb bound states below perturba- 
tive threshold, not on the contour. The difference stems from the presence of the continuum 
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states in the perturbative calculation. The yfz singularity from the continuum states can be 
distorted by the far-away poles. Effectively, the contribution to the imaginary part of the 
correlator in (|5|-[7|) is a product of the real part of the Coulomb bound state poles with the 
imaginary part from the continuum states. 



4 Bottomonium test 

In this section we evaluate the right and left sides of the sum-rule ([TJ]) for bottomonium 
states. We begin with the right side. For our purposes we use the approximation 

G{z, a s ) « a s C F ■ Gi^C* -1 / 2 ). (17) 

where G\ is defined in ([/J), with values m b = 4600 MeV and a^fm^CV = 0.28. The state- 
averaged square of the wavefunction at the origin, |^^(0)| 2 , can be determined from the 
leptonic widths of the s-wave states. The relation is 

-M 2 , , 



l«(o)l 2 ' + ' 



16nQ 2 al ri 



where T e+e _ is the decay width into e + e~, M is the mass of the meson (i.e., 2m + zm), and 
a em is the electromagnetic coupling. 

Let us define E = zm. The left side of ( |T6| ) can be obtained in the following manner. We 
first define a discrete function n(E) such that n(Ei s ) = 1, n(E 2s ) = 2, • • •, where Ek s is the 
energy of the k th s-wave meson. We smoothly interpolate between the values E}. s , thereby 
extending n(E) to a different iable function. We then differentiate n(E) to get the density 
of states, which we call p\{E). We will be working with several versions of the densities of 
states and it is necessary to give them different names. Let us define p 2 {E) as the right side 
of 

m 2 G(z, a s 
4?r 2 1^(0)1' 

In Figure 2 we have plotted pi(E) and p 2 {E) as functions of the meson mass 2m + E. The 
results indicate that the sum-rule is satisfied rather well. 



P^ E ) = , ,, , a „ MM2 ■ ( 19 ) 



5 Determination of the potential 

In this section we propose several new methods for calculating the bottomonium confining 



potential. The sum-rule (16) gave us two calculations for the density of s-wave states, p\ 
and p2- In this section we use WKB methods to convert p± and p 2 into their corresponding 
potentials. 
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Let us write the potential for the heavy quarkonium system as 

V f r ) = -^L + UM. (20) 
r 

We now take a moment to discuss the expected behavior of U(r). In section 4 we found 
that the number of s-wave states with energy less than zm grows more rapidly than any 
inverse power of a s as m — > oo. The number of nodes of the s-wave wavefunctions with 
energy zm must then scale faster than any inverse power of a s . This implies that U(r) is 
slowly varying on distance scales of size m~ 1 z~^, the approximate distance between nodes. 
We also expect that the magnitude of U(r) is small near r = 0. This is because the leading 
term in the perturbative calculation, a s Cp • Gi(a s CFZ~ 1 ^ 2 ), is produced by the continuous 
spectrum of the Coulomb potential. The bound states of V(r) must therefore approximate 
the continuous states of the Coulomb potential and the contribution from U(r) near r = 
must be small. 

For the purposes of calculating energy levels and wavefunctions, we can ignore the change 
of U(r) on length scales Ar < mT 1 and approximate V{r) by Vj(r): 

V I (r) = -?^ + U(r + m- 1 ). (21) 
r 

This approximation introduces an error of size -^Tm^:- 

Since U{m~ 1 ) <C E <C m, we can subtract a constant term, [/(m" 1 ), from V by a small 
redefinition of the heavy quark mass 

m -> m+ ^[/(m -1 ) (22) 

Vj V n = Vi- U{m- 1 ) ■ (23) 

This shift introduces an error on the wavefunctions and masses of the meson states of size 
^Uim^ 1 ). Combining (|23|) with (pl|), we have 



Vj/(r) = -— +W(r), (24) 
r 

where 

W(r) = U(r + 7TI- 1 ) -Uim- 1 )- (25) 

We now apply the WKB approximation. In the WKB approximation the density of states 
is given by 

p(E) = ^— / 26 
* J o 2JE - V n {y) 
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Let us define a function F(x) implicitly as 

1 dW(r) 



F(W(r)) dr ' 
where we have assumed that W(r) is an increasing function of r. We then have 



where the size of the error is n(E) 1 . The Coulomb term effects the integrand of ( p6|) only in 
the small region < y < E^ 1 . As a result, the Coulomb term makes only a small contribution 
to p(E) of size a s . We ignore this contribution and approximate p(E) as 

p(E) » / . ■ 27 

1 ' 2tt Jo JE-W{r) 



(28) 



, . f E F(Y)dY , . 

which we recognize as an Abel transform of F. We can invert the transform, 

v ; ^dEJo VE-Y V ; 

Combining (p8|) and (|30|) and integrating, we find 

2 /-w p(Y)dY , . 

r W = -= / P ; ; 31 

This gives us the inverse of W(r). Using pi and p2 from the previous section, we can calculate 
the corresponding potentials W\ and W 2 - Doing this calculation, we find that W\ and W 2 
agree quite well. The results are shown in Figure 3. Purely for comparison, we plot a s Cp/r, 
the negative of the Coulomb potential, to emphasize that it is not important compared to 
the confining potential in most of the region of interest. 



6 Discussion 



The large- ^ sum-rule provides a new perspective on the heavy quarkonium confining po- 
tential by drawing an explicit relation between the density of states and the wavefunction at 
the origin. As we have seen, the confining potential can now be determined in two different 
ways. In this section we discuss our approximations and consider the question of how to use 
this new information to determine the confining potential more accurately. 

We have made two approximations in the derivation of our sum-rule, large N c and small 
a s . We do not believe that the large iV c approximation introduces a very serious error in this 
analysis. It is used to argue that the bound-state poles are close to the real axis. Because 



the states we actually use in our analysis are rather narrow, this is probably reasonable. By 
working to first order in ot s Cp in ([]), we are ignoring a number of potentially important 
effects, for example the running of the coupling. Thus, for example, we cannot say whether 
we should evaluate the coupling at m 2 , as we have done, or at a smaller mass scale like 
E 2 = (zm) 2 . We are thus ignoring terms proportional to powers of a s ln(l/z). We must do 
this for our perturbation series to be consistent. Because we work to all orders in a s j ^/z, we 
can go down to z for which this is of order one, but for a s ln(l/z) to be large requires still 
smaller z ~ e~ B ^ as ~ AgcWm, which implies E = zm « Aqcd- We expect perturbation 
theory to break down in this region. 

Our calculation of W\ and W2 was done by means of integral transform methods, rather 
than parameter fitting. We believe this carries certain advantages, such as calculational 
simplicity and the elimination of initial biases for the behavior of the potential. Our results 
for W\ and W% appear to indicate a logarithmic or small fractional power law behavior for the 
confining potential in the relevant region for the observed bottomonium spectrum (r < 10~ 2 
MeV -1 ). There are several ways to use our results to produce a more accurate representation 
of the confining potential. One approach is to use W\ or W2 or some average of the two as an 
initial result and then use numerical fitting to find the best possible correction. We expect 
this correction to be quite small, however, since as we now show the initial result appears to 
be quite accurate already. 

We incorporate data for both the s-wave energy levels and leptonic widths in the most 
simple way by constructing an average density of states, 

pUE) = em+m. ( 32) 

We then use ([H]) to find the corresponding potential for p ave (E). As before we have taken 
m& = 4600 MeVand a s (m 2 )CF = 0.28. We now put the Coulomb term back in and calculate 
the resulting spectrum and leptonic widths. The results are shown in Table 1. 



7 Summary 

In our analysis of the heavy quarkonium confining potential we have answered the three 
questions posed in the introduction, and the answer to each question is positive. We derived 
a sum-rule that relates the imaginary part of the heavy quark current-current correlator with 
the product of the density of s-wave states and the state-averaged square of the wavefunction 
at the origin. As we have seen, application of the sum-rule provides a useful complement to 
non-relativistic potential methods. 
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Table 1 



Level 


Mass (th.) 


Mass (ex.) 


r e+e - (th.) 


r e+e - (ex. 


Is 


9490 MeV 


9460 MeV 


1.3 keV 


1.3 keV 


2s 


10040 


10020 


.63 


.52 


3 s 


10370 


10360 


.44 


47 


4s 


10600 


10580 


.34 


.59 


5s 


10820 


10860 


.34 


.31 


6s 


11000 


11020 


.17 


.13 


2p 


9900 


9890 






3p 


10260 


10260 







Figure Captions 

Figure 1. Plot of the integration contour C . 
Figure 2. Plot of the densities p\ and p^. 

Figure 3. Plot of the confining potentials W\ and Wi- The dotted line is a s 
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Figure 1 




Figure 1. Plot of the integration contour C. 
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Figure 2 
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Figure 3 
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Figure 3. Plot of the confining potentials W\ and W 2 . The dotted line is a s Cp/r. 
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